Muon capture on Chlorine-35. 



S. Arole ^, D.S. Armstrong ^, T.P. Gorringe ^, M.D. Hasinoff^, M.A. Kovash^, V. Kuzmin 
^ B.A. Moftah ^t, R. Sedlar T.J. Stocki and T. Tetereva ^ 
^ Department of Physics and Astronomy, University of Kentucky, 
Lexington, KY, USA 40506. 
^ Department of Physics, College of William and Mary, 
Williamsburg, VA, USA 23187. 
^ Department of Physics and Astronomy, University of British Columbia, 
Vancouver, BC, Canada V6T IZl. 
^ Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, Dubna, I4I 

980, Russia. 

^ Dubna Branch of Skobeltsyn Institute of Nuclear Physics, Lomonosov Moscow State University, 

Dubna, I4I 980, Russia. 
(February 8, 2008) 



Abstract 

We report measurements of 7-ray spectra from muon capture on ^^Cl. For 
the allowed Gamow-Teller transitions to the ^^8(2939,3/2"*") state and the 
^^8(3421,5/2"*") state we obtained their capture rates, hyperfine dependences 
and correlation coefficients. The capture rates and hyperfine dependences 
were compared to shell model calculations using the complete ls~Od space 
and the universal 8D interaction. With gp/ga = 6.7 and ga = —1.00 (or 
ga = —1.26) we found agreement of the model and the data at the 1~2 a 
level. However, we caution that the transitions are sensitive to i = 2 forbidden 
matrix elements. 
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I. INTRODUCTION 



The induced pseudoscalar coupling (gp) is the least known of the proton's 
weak couplings. For the free proton, the coupling's determination is an im- 
portant test of chiral symmetry breaking For the bound proton, the 
coupling's renormalization is sensitive to vr exchange currents, A-hole excita- 
tions and possible precursor effects of chiral phase transitions in hot, dense 
nuclear matter [Q-^. 

Unfortunately the effects of gp are subtle and elusive. In few-body systems 
recent results from radiative muon capture (RMC) on and ordinary 

muon capture (CMC) on ^He are available. However their interpretations 
are complicated by muon chemistry in and 2-body currents in ^He, and the 
puzzling discrepancy between the values of gp from RMC on and OMC 
on ^He is so far unexplained. In complex nuclei new results from allowed 
transitions on ^^B [^], ^'^Na [11| and ^^Si p!^-p!^ are also available. In nuclei 
the difficulty is disentangling the weak dynamics from nuclear structure. The 
majority of data on nuclei are consistent with an unrenormalized gp, but more 
experiments on other transitions would be interesting. 

In this article we report measurements of 7-ray spectra from muon cap- 
ture on ^^Cl. In particular we describe the determination of capture rates, 
hyperfine dependences and 7-1/ angular correlations for two ^^Cl — > ^^S al- 
lowed Gamow-Teller transitions. We also compare their capture rates and 
hyperfine dependences to a large basis shell model calculation, and discuss 
their sensitivity to the weak couplings and the nuclear structure. 

The paper is organized as follows: In sections || and ID we describe the 
measurement setup and experimental results. In section IV we discuss the 
comparison of the model and the data, and the sensitivity to the couplings 
constants and the nuclear model. We conclude in section 



II. EXPERIMENTAL SETUP 

The experiment was conducted using the M9B backward decay muon 
beamline at the TRIUMF cyclotron. The setup was similar to our earlier 
studies of other Is-Od shell nuclei (see Refs. 111,13,151). 



We employed a /u~ beam of incident momentum 60 MeV/c, stopping rate 
1.2 X 10^ s~^, and e and vr contamination of 10% and <0.1%. Muon stops 
were counted in a plastic scintillator beam telescope comprising two counters 
(SI and S2) upstream of the target and one counter (53) downstream of the 
target. The target material was isotopically enriched Na^^Cl (99% chlorine- 
35) powder of mass 50 g (we used isotopically pure chlorine-35 to avoid 7- 
ray backgrounds from '^'^Cl{fi,2ni^) reactions). The material was packed in a 
thin-walled, disk-shaped, polyethylene container. The target was viewed by 
a HPGe detector with an in-beam energy resolution of 2.2 keV (full width 
half maximum) and in-beam time resolution of 10 ns (FWHM) for 1.33 MeV 
^'^Co 7-rays. A segmented Nal annulus, surrounding the HPGe detector, was 
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used to suppress the Compton scattering background. A second Nal array, 
viewing the Na^^Cl target, was used for cohection of 7-ray coincidence data. 

Events were digitized on fulfillment of the logic condition Ge-CS-^STOP 
■busy, where Ge indicates a signal in the HPGe detector, CS indicates no 
signal in the Nal suppressor, ^STOP indicates a fi^ stop in the preceeding 
2.0 /is, and busy indicates that the acquisition is live. For each event we 
recorded energy and timing signals from the Ge detector, Nal arrays and 
beam counters. We also recorded a multi-hit time history of muon stops, and 
a pile-up bit and an overload bit from the HPGe. The pile-up bit identifies 
events with a preceeding Ge detector hit within 50 //s and the overload bit 
identifies events with a preceeding > 10 MeV Ge detector hit within 500 /is 
(for details see Ref. [^). Data were collected for 4.2 x 10^'^ muon stops in 
the Na^^Cl target, along with 7-ray background data from a LiCl target and 
x-ray calibration data from various materials. 

Offline the data were sorted into histograms corresponding to (i) a HPGe 
singles 7-ray energy spectra, (ii) a HPGe-Nal coincidence 7-ray energy spec- 
tra (i.e. a HPGe energy spectrum with a Nal coincidence requirement), and 
(iii) many MPGq-^STOP time-binned 7-ray energy spectra {i.e. HPGe en- 
ergy spectra corresponding to different time bins). The time bins were defined 
by the time difference between the incoming muon and the outgoing gamma- 
ray. In filling the histograms we rejected any events with either an overload 
signal or a pile-up signal in the HPGe. 



III. EXPERIMENTAL RESULTS 
A. HPGe calibration 

A large body of muonic x-ray data was collected in order to determine the 
HPGe photon acceptance, energy resolution and time resolution. The data 
comprised K, L and M series x-ray spectra from S, Ca, In, Nb, Fe, Pb and 
Bi, and spanned a range of energies from 500 to 4000 keV. 

To determine the acceptance with muonic x-rays we either (i) used pub- 
lished yield data for individual x-rays ||T6|"[l^ or (ii) assumed a yield of unity 
for the entire K-series. Then we performed a least squares fit of the mea- 
sured acceptances to a smooth empirical curve to permit interpolation to the 
interesting 7-ray energies. 

For the energy resolution determination we fit the numerous x-ray peaks 
to a central Gaussian supplemented by both low and high energy tails. The 
width of the central Gaussian was observed to increase from 2.2 keV at 
1000 keV to 3.6 keV at 3000 keV. Both low energy and high energy tails 
were observed, but their effects on the fits to the ^^S 7-rays were almost neg- 
ligible (see Sec. HI G| ). The best fit values of the lineshape parameters were 



used to fix the instrumental lineshape in the least squares fits to the 7-ray 
data. 
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For the time resolution determination we fit the prompt peaks of the x- 
ray hnes with a central Gaussian supplemented by both early and late tails 
(energy windows were used to select the x~rays). We found that, for photon 
energies above 1000 keV, the Gaussian width was approximately constant at 
10 ns FWHM, but the resolution rapidly got worse for lower energies. A weak 
tail at late times was observed, but its effect on the fits to the ^^S 7-rays was 
almost negligible (see Sec. Ill F| ) . Again, the best fit values of the lineshape 
parameters were used to fix the lineshape in the least squares fits to the 7-ray 
data. 

More details on the procedures for the determination of the photon ac- 
ceptance, energy resolution and time resolution can be found in the Refs. 

niPl- 



B. Line Identification 

About eighty peaks were found in the 7-ray spectrum from muon stops 
in Na^^Cl. To match the peaks to known 7-ray lines we searched the nu- 



clear 7-ray database at the National Nuclear Data Center (NNDC) [20| for 
all isotopes corresponding to the capture reactions ^^Na(/U~ , xp yn u) and 
^^Cl(//~,xp yn u) where x,y = 0,1,2. For conclusive identification we de- 
manded (i) an energy match within experimental uncertainties, (ii) consis- 
tency with other 7-branches of the parent state, and (iii) Doppler broadening 
if the parent state lifetime was less than the recoil stopping time. This pro- 
cedure yielded 5 7-rays from levels in '^^S, 21 7-rays from levels in ^^S, and 
4 7-rays from levels in '^^S. 

The conclusively identified 7-rays from ^^S are listed in Table | and shown 
m Fig. 1. The five 7-rays all feed the ^^S ground state and correspond to three 
allowed transitions to the (2717,5/2+), (2939,3/2+) and (3421,5/2+) states 
and two forbidden transitions to the (1991,7/2") and (2439, 3/2~) states. 

In addition, the production of 7-rays from the '^^8(1572,1/2+) and 
35S(4028,l/2+-5/2+) levels is possible. For the (1572,1/2+) level, a 7-peak 
was clearly identified at 1572 keV, but unfortunately it matches the energies 
of both the 1572 ^^S transition and a 4877 3304 transition. The 
observed 1572 keV peak is Doppler broadened, indicating that the short life- 
time ^^S state is the major 7-ray source. However some contribution from 
the ^^8(1572, 1/2+) level is a possibility. For the (4028, 1/2+-5/2+) level, 
weak lines were identified at energies of -^5 keV below that expected for 
4028 ^ 1572, 4028 2348 and 4028 2939 transitions. Note the experi- 
mental uncertainty for the 4028 keV excitation energy is about ±2 keV p§. 



so these peaks in the Na^^Cl data may be evidence that the energy of the 
state is actually about 4023 keV. 
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C. Line intensities 



To determine the gamma-ray yields per fiCl atom formed, we employed 
the equation 

= eAniE) ^'^^ ^'"'^ 

where is the number of 7-rays detected, is the number of livetime- 
corrected muon stops, eAQ,{E) is the photon acceptance at the appropriate 
energy, fiso is the fi~ CI atomic capture fraction, and Cab and Csv are minor 
correction factors, described below. The 7-ray counts were obtained from 
fits with Doppler broadened or Gaussian lineshapes (see section III A| ). The 



atomic capture fraction fiso = 0.59 =b 0.04, i.e. the fraction of muon stops 
in NaCl that undergo atomic capture on CI, was taken from Ref. The 
factors Cab and Cgv accounted for photon absorption in the target and self- 
vetoing by the suppressor {Cab varied from 1.02 to 1.03 and Csv varied from 
1.06 to 1.16). They are discussed in detail in Ref. [pT|| . 

The resulting gamma-ray yields per /iCl atom formed are given in Table |. 
The dominant uncertainties were the measurement uncertainty in the atomic 



capture fraction of ±7% (see Ref. ||21[] for details) and the normalization un- 



certainty in the muonic x-ray calibration of ±15% (see Ref. ||TT| for details). 
The statistical uncertainties in and total uncertainties in Cab and Csv were 
negligible. 



D. Cascade feeding 

A serious concern in 7-ray studies is cascade feeding into interesting ^^S 
levels from higher-lying levels. If missed, such production of ^^S states would 
distort the interpretation of the measured rates, hyperfine dependences and 
angular correlations. 

Although the comparison between our Na^^Cl data and the NNDC Tables 
|20| gave only five clean 7-ray matches, some 7-peaks were never identified 



and many ^^S states have unknown decays. Therefore we performed a 7- 
7 coincidence measurement in order to study the total amount of cascade 
feeding to ^^S levels. Specifically, we determined the counts for interesting 
7-ray lines in the HPGe singles spectrum (denoted N^) and the HPGe-Nal 
coincidence spectrum (denoted Nj), and computed the super-ratio 



f- ^y^'s 



For comparison the Fermi Z-law yields fis, 
fiso = 0.58. The earlier experimental work of Knight et al. |22| gave fiso 
gave fiso = 0.49. 



0.61 and the calculation of Vogel et al. yields 

0.56 and Zinov et al. 
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The coicidence fraction fc reflects the total number of other 7-rays in prompt 
coincidence with the interesting 7-ray. Note in Eqn. ^ the ratio Nq°^^ /Ng°^^ 
was obtained by measuring the well known 1.17 and 1.33 MeV coincident 
7-rays from a Co-60 source. The ratio Nq°^^ /N^"^^ serves to normalize the 
ratio N^/Ng. We stress that the approach has limitations: (i) it cannot 
distinguish feeding via a single-step cascade from a multi-step cascade, and 
(ii) its interpretation assumes an energy-independent acceptance for the 7- 
ray detection by the Nal array (the latter assumption is fairly reasonable for 
7-ray energies from 0.5 to 3.0 MeV). The method is described in detail in 
Refs. lll|,|l|. 

The measured coincidence fractions fc for ^^S 7-rays are listed in column 
seven of Table |. For the 1991, 2348 and 2717 keV 7-rays the table shows the 
presence of considerable feeding from unidentified levels. The comparison of 
experiment and theory is therefore not warranted in such cases. However for 
the 2939 and 3421 keV 7-rays the table shows the absence of large amounts 
of cascade feeding. Therefore we focused our attention on the model-data 
comparison of the physical observables in the ^^01(3/2+, 0) 253(3/2+^ 2939) 
transition and the 35Cl(3/2+,0) ^^8(5/2+, 3421) transition. 



E. Capture rates 

To convert the gamma yields to capture rates we accounted for the 7-ray 
branching ratios (taken from Rcf. ||2^) and multiplied by the muon disap- 
pearance rate (taken from Ref. For 35ci(3/2+,0) ^ 253(3/2+, 2939) 
and 2^01(3/2"'", 0) ^^8(5/2+, 3421) the resulting rates and experimental un- 
certainties are listed in Table 

Because of the hyperfine effect in the fi~Cl atom, the experimental capture 
rates are combinations of hyperfine capture rates (for details see 8ec. IIIF). 



Assuming an initial statistical population of the two hyperfine states, the 
hyperfine transition rate of 8ec. [II F| , the published muon disappearance rate 



of Ref. and a 2.0/US wide //^8T0P gate, the observed rates A in Table 
^ are related to the hyperfine capture rates via A = 0.12A+ -|- 0.88A_. The 
extraction of the hyperfine capture ratio A-(_ /A_ is discussed in the following 
section. 



F. Hyperfine dependences 

In a non-zero spin / 7^ target the muonic atom's 18 atomic state is 
split into hyperfine states with F_|_ = / -|- 1/2 and F_ = I — 1/2. Due to the 
spin dependence of the weak interaction the capture rates from the hyperfine 
states are usually different. In certain cases {e.g. ^^Cl), during the muonic 
atom lifetime the upper HF state decays to the lower HF state. Consequently, 
the time dependence of 7-rays from ^~ capture has the form 
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where A/) is the disappearance rate, Ah is the hyperfine transition rate, and 
k is related to the hyperfine dependence of muon capture via k = (A-|- /A_ — 
1). The quantity /+ is the fractional population of the state at t = 0. For 
a spin 3/2 target and a statistical HF population |26| the factor is /+ = 
5/8. Consequently, if Ah is neither too fast nor too slow {i.e. Ah ^ Ad), 
the hyperfine dependence of muon capture can be determined from the time 
spectrum of the 7-rays. 

For each 7-ray line its time-binned energy spectra were first fit to appropri- 
ate Gaussian or Doppler lineshapes in order to determine its time spectrum. 
The 7-ray time spectra themselves were then fit to a convolution of the the- 
oretical time dependence (Eqn. with the measured HPGe time resolution 
(Sec. Ill A| ). We accounted for the slight distortion of the disappearance rate 
due to muon pile-up in the Na^^Cl target by the method described in Ref. 
Ijllj (a 4% effect). Additionally we investigated the sensitivities of A+/A_ 
and A/j to the various parameters of the instrumental response and the muon 
pile- up. 

In this analysis we employed the 1991 keV 7-ray time spectrum to deter- 
mine the hyperfine transition rate Ah (since the 1991 keV 7-ray has a very 
large hyperfine effect). Note, our best fit value of A/j = lO.litl.O /is~^ for this 
NaCl experiment and the earlier value of A/j = 8.1 ±2.2 fj.s~^ for a LiCl exper- 
iment 1 27 1, are consistent (we are also consistent with the recent experiment 
of Stocki et al. [|2^). Also our experimental rate is in reasonable agreement 
with the predicted rate of Ah = 8.0 /us~^ from Ref. |29|. 

The best fit values for the hyperfine dependences of the 2939 keV line and 
3421 keV line are given in Table m (in extracting A+/A_ we fixed A/j using the 
1991 keV 7-ray). In addition, representative fits to their 7-ray time spectra 
are shown in Fig. |2| and representative sensitivities to the various fitting 
parameters are shown in Tables III and |^ The quoted errors for A_|_/A_ 
include the statistical errors in the fitting procedure and the sensitivities to 
the determination of the instrumental resolution and the muon pile-up. 



G. Angular correlations 

Because of the spin dependence of the weak interaction, the recoil pro- 
duced in muon capture is usually oriented. In general, this recoil orientation 
yields an angular correlation between the neutrino direction and the decay 
7-ray direction. Further, when the 7-decay lifetime is less than the recoil 
stopping time, the ^-v directional correlation will be manifest in the 7-ray 
Doppler spectrum. Both the 2939 keV 7-ray and the 3421 keV 7-ray were 
Doppler broadened. 

For unpolarized muons the directional correlation of the photon and the 
neutrino has the form 



1 + aP2{cos9), 



where 9 is the angle between the momentum vectors of the photon and the 
neutrino, P2{cos9) is the Legendre polynomial, and a is the ^-v angular 
correlation coefficient. Note, for recoil spins J < 2 or 7-ray multipolarities 
L < 2, higher-order Legendre polynomials do not contribute. 

To fit the 7-ray Doppler spectra, and extract the ^-v correlation coef- 
ficient, we convoluted the theoretical Doppler spectrum with the HPGe in- 
strumental lineshape. In the fitting procedure we fixed the lineshape param- 
eters at their known values, and studied the background sensitivity by trying 
background shapes of various forms {e.g. linear, quadratic and exponential). 
Representative fits to the 2939 keV and 3421 keV lines are shown in Figs. ^ 
and ^ and typical sensitivities to the fitting parameters are listed in Tables 
and For the 2939 keV line a concern is clearly the background line at 



2935 keV, which obscures a portion of the Doppler broadened spectrum of 
the 2939 keV gamma-ray. 

Unfortunately the Doppler lineshape may be distorted by the slowing 
down of the recoil ion. Specifically, if the slowing-down time tg and the 
gamma-ray lifetime r are comparable, then the lineshape is a function of 
the coefficient a and the ratio r/ts- From sample fits with different lifetimes 
we found the input value of r/ts and output value of a were highly correlated 
when r/ts > 0.05 (see Figs. | and |). For r/ts < 0.01 the 2939 keV line 
yielded a = -0.43 ± 0.13 and the 2939 keV line yielded a = -0.39 ± 0.05 
(the quoted errors for the coefficients include the statistical uncertainties from 
the fitting procedure and the various sensitivities to the fitting parameters). 
However, for r/ts > 0.01 the magnitude of a decreases as the value of r/tg 
increases. Therefore, since independent determination of the gamma-ray life- 
times are currently unavailable, our results for a are functions of r/tg {i.e. 
the plots of Fig. ^ and Fig. ^). When the necessary lifetimes are finally 
measured the correlation coefficients may be read off these figures. 

Finally because of the hyperfine effect, the observed coefficient a is a linear 
combinations of the two coefficients a± of the two F± states. Specifically the 



observed correlation a is (see Ref. |37| for details) 



n+A+a+ -|- n_A_a_ 

= 1 1 ) 

n+A-i- -|- n_A„ 

where n^/n^ is the relative muon occupancy of the hyperfine states and 
A-|-/A_ is the relative capture rates from the hyperfine states. The relative 



occupancy was taken from Sec. IIIF, and is n+/n_ = 0.12/0.88 = 0.14 



The hyperfine dependences were taken from Table |l|, and are 0.74 it 0.17 for 
the (2939,3/2+) transition and 1.60 ± 0.19 for the (3421,5/2+) transition. 
This yields a = 0.09a+ + 0.91a_ for the (2939,3/2+) transition and a = 
0.18a+ + 0.82a- for the (3421,5/2+) transition. 
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IV. INTERPRETATION 



Herein we discuss the comparison of experiment and theory for the mea- 
sured observables in the ahowed transitions to the (2939,3/2+) level and the 
(3421, 5/2"'") level. In Sec. IV A we describe the general dependences on weak 
dynamics, and in Sec. [VB| we discuss the model calculation and input param- 
eters. Sec. IV C compares the model results with the experimental data, and 
emphasizes the sensitivities to the weak couplings and the nuclear structure. 

Note that we restrict the comparison of theory and experiment to the cap- 
ture rates and the hyperfine dependences. Interpretation of the correlations 
is not possible at the moment as the 7-decay lifetimes and their mixing ratios 
are not currently available. 



A. Fujii-PrimakofF Approximation 

To demonstrate the basic features of the physical observables in allowed 
3/2"*" 3/2+ and 3/2+ — > 5/2"*" transitions we first review some results 



of the Fujii-Primakoff approximation [31|. Recall that the Fujii-Primakoff 



Hamiltonian for nuclear muon capture [32| is 
I - (J ■]} ^ 

H = r+ ^ T-"(Gi/ I -li + GA CF ■ Gi + Gp a ■ v at - u)5{r - r^), 

^ i=i 

where v is the z/-momentum unit vector, 1 and a are unit and spin matrices 
(the operators with subscripts act on nucleons and the operators without sub- 
scripts act on leptons), and t+ converts the muon into a neutrino. Note that 
the leading contribution of ga is to the effective coupling Ga and the lead- 
ing contribution of Qp is to the effective coupling Gp. In the Fujii-Primakoff 
approximation only the allowed GT operator and the allowed Fermi operator 
are retained following the multipole expansion of Eqn. |^. 

First consider the hyperfine dependence of muon capture. We note in the 
multipole expansion of the Fujii-Primakoff Hamiltonian the G^-term makes 
allowed contributions to neutrino waves with total angular momentum j'^ = 
1/2+ only whereas the Gp-term makes allowed contributions to neutrino 
waves with total angular momentum j'^ = 3/2+ (this difference is because of 
z> in Eqn. In 3/2+ — > 5/2+ transitions, neutrinos with j = 1/2+ may be 
emitted in capture, but neutrinos with j = 3/2+ must be emitted in F_ 
capture. This makes A-|-/A_ >> 1 and a strong function of the ratio Qp/ ga- 
However in 3/2+ 3/2+ transitions the emission of j = 1/2+ neutrinos is 
possible for both atoms and F- atoms. This makes A_|_/A_ ~ 1 and a 
weaker function of the ratio Qp/ Qa- 

Next consider the 7-2/ correlation in muon capture. Note in Eqn. |^ the 
Gyi-term is a vector in spin-space while the Gp-term is a scalar in spin-space. 
Consequently, in most cases the magnetic sub-states of recoil nuclei are pop- 
ulated differently by allowed contributions that originate from the G^-term 
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and the Gp-term. This makes the recoil orientation, and therefore the ■j-v 
correlation, a function of Qp/ga- However an exception is ^-v correlations in 
3/2"'" — > 5/2"*" transitions on F- atoms. In such circumstances the emission 
of d-wave i/'s is required and only the Gp-term generates an allowed contri- 
bution. Consequently for 3/2"*" ^5/2"*" transitions on F_ atoms the recoil 
orientation and ^-v correlation is independent of the coupling constants in the 
Fujii-Primakoff approximation. The potential usefulness of the 7-1/ correlation 
coefficients for model testing is illustrated further in Appendix 

In summary, different observables offer different sensitivities to the weak 
couplings and the nuclear structure. For example, the hyperfine dependence 
A-|-/A_ in a 3/2"*" ^5/2^ transition is particularly sensitive to the induced 
pseudoscalar coupling and the correlation coefficient a- in a 3/2"*" —5-5/2" 
transition is particularly sensitive to the forbidden matrix elements. In prin- 
ciple they enable a means of both extracting the weak couplings and testing 
the nuclear structure. 



B. Model calculation 



Our model calculations of the capture rates and hyperfine dependences 
were performed in the impulse approximation using the shell model (equa- 
tions for the capture rates and the hyperfine dependences are published in 
Refs. l3|,||,|§). Specifically, we used the computer code OXBASH |3|, the 



complete Is-Od space, and the universal SD interaction ^6|. We fixed the 
weak vector and magnetic couplings to the values = 1.000 and gm = 3.706 
and varied the weak axial and induced pseudoscalar couplings. The j4 = 35 
nuclear matrix elements were computed with harmonic oscillator wavefunc- 
tions and an oscillator parameter h = 1.90 fm. The momentum transfer 
was computed using the experimental values of excitation energies (giving 
q = 101.17 MeV/c for the 2939 keV state and q = 100.59 MeV/c for the 
3421 keV state). Finally the /U~^^C1 atomic wavefunction was assumed to be 
uniform in the nuclear volume and computed employing a muon wavefunction 
reduction factor of i? = 0.521 (see Walecka |3^ and references therein). More 
details are given in Ref. |jll|. 



The arguments we made in Sec. IV A for sensitivities of observables to 



ga and gp were based on the dominance of the allowed GT matrix element 
in the interesting [ji, v) transition. However our model calculations show a 
large (0^5/2)^^(181/2)^(0^3/2)^ component in the ^^Cl ground state and a large 
(0^5/2)^^(151/2)^(0^3/2)'* component in the two ^^S states. Consequently the 
1^1/2 ~^ 0^^3/2 single particle transition is important in both ^^Cl(3/2+,0) 
^ 35s(3/2+,2939) and 35ci(3/2"^, 0) ^ ^^S{b/2+ ,U2l). Furthermore for 
1^1/2 ~^ 0^^3/2 single particle transitions the allowed GT matrix element is 
zero (the operator obeying a A£ = selection rule). This amplifies the im- 
portance of £ = 2 forbidden contributions in both transitions. A serious 
concern is clearly therefore the correct accounting for 1 = 2 forbidden contri- 
butions from the lsi/2 single particle transition in the ^^Cl(3/2+,0) 
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^ 35s(3/2+,2939) and the 3^01(3/2+, 0) ^ 353(5/2+^3421) |]. 

We further note that the calculation shows that the ^^Cl ^ ^^S GT 
strength is spread over many states with Ex < 10 MeV (see Fig. ^). This 
contrasts with ^^Na — > ^^Ne and ^®Si — >■ ^®A1 (other cases of /i capture work 
on Is-Od nuclei) where only a few states are found to exhaust a large fraction 
of GT strength (for more details see Gorringe et al. [p!5| ). 



C. Comparison of experiment and theory 

Fig. 1^ shows the results of our calculations of the capture rate A = 
0.12A+ + 0.88A_ to the (2939,3/2+) state and the (3421,5/2+) state. As 
expected we found that the rate is most sensitive to the value of the coupling 
Qa- For ga = —1.26 and gp/ga = 6.8 the computed rates were 12.2 x 10'^ s~^ 
for the (2939,3/2+) state and 15.4 x 10^ for the (3421,5/2+) state. For 
Qa = —1.00 and Qp/ga = 6.8 the computed rates were 8.4 x 10'^ s~^ for 
the (2939,3/2+) state and 12.3 x 10^ s"^ for the (3421,5/2+) state. The 
calculated rates are remarkably close to the corresponding measurements of 
(12.2 ± 2.2) X 10^ s-i and (11.9 ± 2.2) x 10^ s"^ respectively 

Concerning the induced pseudoscalar coupling, the calculated (3421, 5/2+) 
rate was found to exhibit some sensitivity to gp but the calculated (2939, 3/2+) 
rate was found to exhibit no sensitivity to gp (see Fig. ^). This finding is in 
agreement with the expectations of the Fujii-Primakoff Approximation (see 



Sec. [VA) 



Figs. ^ show the results of our calculations for the hyperfine dependence 
A+/A_ of the (2939,3/2+) transition and the (3421,5/2+) transition. Note 
as argued in Sec. p^V A| , the 3/2+ — > 5/2+ hyperfine effect is quite strongly 
dependent on gp whereas the 3/2+ — > 3/2+ hyperfine effect is quite weakly 
dependent on gp. For gp/ga = 6.7 the calculation gives A_|_/A_ = 1.3-1.4 
for the (3421,5/2+) transition and A+/A_ = 0.45-0.46 for the (2939,3/2+) 
transition (the range corresponds to choosing either ga = —1.00 or ga = 
— 1.26). For comparison the experimental values are 1.60±0.19 and 0.74ib0.17 
respectively. 

Clearly the model calculations and experimental results for hyperfine de- 
pendences are in semi-quantitative agreement; for example both concurring 
that -F_ capture is stronger for the (2939,3/2+) transition and -F+ capture 
is stronger for the (3421,5/2+) transition. However, at the level of 1-2 o", 
some indication of disagreement between model and data is suggested for 
Qp/ Qa = 6.7, the data favoring a smaller value of gp/ga- 



^The full calculation shows numerous other configurations with additional holes in the 0^5/2-153/2 
orbitals of the relevant ^ = 35 states. These generate single particle contributions in 3^Cl(3/2+,0) 
3^S(3/2+, 2939) and ^^01(3/2+, 0) 3^S(5/2+, 3421) involving Od ^ Off and Is Is transitions. 
Such single particle transitions generate substantial allowed Gamow- Teller contributions. 
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Is this evidence for the medium modification of the induced couphng? 
Recall in discussing the A = 35 structure in Sec. IV B we identified the 
large contribution of lsi/2 — > 0^3/2 transitions and i = 2 forbidden operators. 
For example, the fact that A+ ~ A_ not A+ >> A_ for ^^Cl(3/2+ , 0) 
^^S( 5/2+ , 3421), both in the data and in the model, is supporting evidence 
that forbidden terms are important here. This means that the observables 
are sensitive to the ratio between the allowed Gamow-Teller matrix elements 
and the various £ = 2 forbidden matrix elements. The matrix elements are 
rather dependent on the particular admixtures of the various hole-states in 
the leading (0(i5/2)^^(lsi/2)^(0'^3/2)'^ configuration for ^^Cl and the leading 
(0(i5/2)^^(lsi/2)^(0'^3/2)^ configuration for ^^S. Therefore we believe it is un- 
wise to blame the coupling gp for the small discrepancies between the model 
and the data. 

We should also remind the reader that some cascade feeding into the 
(2939,3/2+) state and the (3421,5/2+) state is possible (see Sec. IHTdD . If 
feeding is present, the measured capture rate will over-estimate the direct 
capture rate, and the measured hyperfine dependence may differ from the 
direct hyperfine dependence. This also could account for some remaining 
discrepancies between experiment and theory. 



V. CONCLUSION 

In summary, we have measured the 7-ray spectra from muon capture on 
isotopically enriched Na^^Cl. For the 35ci(3/2+ , 0) ^^S{ 3/2+ , 2939) tran- 
sition and 3^01(3/2+ , 0) ^^S{ 5/2+ , 3421) transition we obtained their 
capture rates, hyperfine dependences and j-v correlation coefficients. Con- 
cerning the capture rates and hyperfine dependences the experimental results 
are in semi-quantitative agreement with a large basis shell model calculation 
using Qp/ga = 6.7 and ga = —1.00 (or ga = —1.26). However we note that 
large £ = 2 forbidden contribution from large lsi/2 ~^ single particle 

transitions are probable indications of model dependences. We therefore are 
unwilling to claim a determination of gp nor ga, with any precision. 

For the angular correlations, since independent measurements of gamma- 
ray lifetimes and mixing ratios are currently unavailable, the comparison of 
model and data was thwarted. With such supplemental 7-ray data the cor- 
relation coefficients for the two 7-rays would permit additional testing of the 
calculation and the extraction of the coupling gp. For further details see Sec. 
IIIG and Appendix ^ We therefore encourage any future efforts to measure 



these quantities. 

We would also encourage a more thorough theoretical investigation of 
these allowed Gamow-Teller transitions. In particular, a more quantitative as- 
sessment of the model uncertainties and the various approximations is worth- 
while. For example, in our model we have employed a uniform muon wave- 
function and harmonic oscillator nuclear wavefunctions, and the effect of these 
simplifications on the observables should be studied. 
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APPENDIX A: USEFULNESS OF THE -f-iy CORRELATIONS. 

In this appendix we describe the usefulness of the j-v correlation coeffi- 
cients in the testing of the nuclear model calculation. We consider the example 
of the 3421 keV gamma-ray from the (5/2-I-, 3421) excited state. Assuming 
negligible slowing-down effects the measured 7-2^ correlation coefficient was 
found to be a = -0.39 ± 0.05. 

In our toy model for the 35ci(3/2+,0) ^ 353(5/2+, 3421) ^ 353(3/2+, 0) 
sequence we will assume that the '^^S{ 5/2"*" , 3421) lifetime is fast enough 
that slowing-down effects are negligible and that the 3421 keV gamma-ray is 
pure Ml radiation. The purpose of the toy model is to illustrate the potential 
sensitivities of the correlation a, and the importance of a measurement of the 
35S(5/2+,3421) lifetime and the 3421 keV E2/M1 mixing ratio. 

According to Ciechanowicz and Oziewicz the 7-1^ angular correlation 
coefficient a may be written as the product 

Q = 02-62 

where 02 is determined by the fi capture process and B2 is determined by 
the 7-decay process. We obtained 02 using the Is-Od shell model with the 
universal SD interaction as discussed in Sec. IIV B|. We took B2 = +0.3741 



from Table 1 of Ciechanowicz and Oziewicz under the assumption of a 
pure Ml decay. Note that we computed the two correlations coefficients 
for the two hyperfine states F± and then combined the values to obtain the 



observed correlation a = 0.18a+ + 0.82a_ (see Sec. pl G| for details). 

Our toy model results are given in Table [VI ]| for ga = —1.26. It shows 
results from the full calculation as well as from a calculation omitting the 



interesting second forbidden matrix elements. As argued in Sec. pV A the 



correlation coefficient is rather weakly dependent on the induced coupling gp 
but somewhat more dependent on second forbidden terms. 

Our measured value a = — 0.35±0.05 and the calculated values a ~ —0.23 
are significantly different. Unfortunately in the absence of any experimental 
data on the state lifetime and the E2/M1 mixing ratio it's impossible to know 
if the discrepancy is a reflection of problems in the model calculation of the 
muon capture or the invalidity of our model assumptions of fast 7-decay and 
pure Ml radiation. We therefore strongly encourage the measurements of r 
and 6 for the ^^S{ 5/2+ , 3421) ^ 35g( 3/3+ ^ g) gamma-decay. 
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TABLES 

TABLE L Gamma-ray yields per /xCl atom formed for the cleanly identified lines from 
^^Cl(//,i^)^^S. Columns 1-5 give the relevant energies, lifetimes and branching ratios from Ref. 
[20]. The coincident fraction fc is discussed in Sec. HID. 





( E, J- ), 




B.R. 


lifetime 


7-ray yield 






coinc. frac. 


(keV) 


(keV) 


(keV) 


(%) 




(per /xCl atom) 




fc 


1991.3 


(1991,7/2") 


(0,3/2+) 


100 


(1.02 ±0.05)ns 


( 7.1±1.2 ) ; 


<10" 


-3 


1.27 ± 0.06 


2347.8 


(2348,3/2-) 


(0,3/2+) 


75 


0.89±0.12fs 


( 2.8±0.5 ) ; 


<10" 


-3 


0.77 ± 0.15 


2717.1 


(2717, 5/2+) 


(0,3/2+) 


94 


(70 ± 25)fs 


( 2.9±0.5 ) ; 


<10" 


-3 


1.00 ± 0.39 


2939.6 


(2939,3/2+) 


(0,3/2+) 


100 




( 5.4±0.9 ) ; 


<10" 


-3 


<0.30 


3421.0 


(3421,5/2+) 


(0,3/2+) 


100 


<70/s 


( 5.3±0.9 ) ; 


<10" 


-3 


0.10±0.15 



A 



TABLE IL The measured capture rates and hyperfine dependences for the '^^Cl(3/2+, 0) 
35S(3/2+,2939) transition and the 35ci(3/2+,0) ^ ^53(5/2+, 3421) transition. Note that 
= 0.12A+ + 0.88A_ (see Sec. for details). 



( E, J- )i 


A 


A+/A- 


(keV) 


(xlO^ s-i) 




(2939,3/2+) 


12.2 ±2.2 


0.74 ±0.17 


(3421,5/2+) 


11.9 ±2.2 


1.60 ±0.19 
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TABLE III. Results of our studies of the systematic uncertainties in the extraction of the 
hypcrfinc dependence for the '^^Cl(3/2+ , 0) '^"''8(3/2+ , 2939) transition. The input parameters 
in the time spectrum fits were the t = position (to), instrumental width (a), muon disappearance 
rate (Ad) and hyperfine transition rate (A/j). See text for details. 



Aft 


Ad 


to 


a 


k 


A+/A- 


xlO^ s-i 


xlO^ s-i 


chans. 


chans. 







8.0 


0.225 


1504.04 


1.923 


-0.15±0.10 


0.76±0.16 


10.0 


0.225 


1504.04 


1.923 


-0.15±0.10 


0.76±0.16 


8.0 


0.235 


1504.04 


1.923 


-0.18zb0.10 


0.71±0.16 


10.0 


0.235 


1504.04 


1.923 


-0.19±0.10 


0.70±0.16 


10.0 


0.225 


1503.84 


1.923 


-0.14±0.10 


0.78±0.16 


10.0 


0.225 


1504.24 


1.923 


-0.16±0.10 


0.74±0.16 


10.0 


0.225 


1504.04 


1.82 


-0.15±0.10 


0.76±0.16 


10.0 


0.225 


1504.04 


2.02 


-0.15±0.10 


0.76±0.16 


10.0 


0.225 


1504.04 


no tail 


-0.15±0.10 


0.76±0.16 



TABLE IV. Results of our studies of the systematic uncertainties in the extraction of the 

hyperfine dependence for the '^^Cl(3/2+ , 0) ^^8(5/2"'' , 3421) transition. The input parameters 
in the time spectrum fits were the t = position (to), instrumental width (a), muon disappearance 
rate (Ad) and hyperfine transition rate (A^). See text for details. 



Ah Ad to a k A+/A_ 

xlO^ s"^ xlO^ s~^ chans. chans. 



8.0 


0.225 


1504.04 


10.0 


0.225 


1504.04 


8.0 


0.235 


1504.04 


10.0 


0.235 


1504.04 


10.0 


0.225 


1503.84 


10.0 


0.225 


1504.24 


10.0 


0.225 


1504.04 


10.0 


0.225 


1504.04 


10.0 


0.225 


1504.04 



1.923 


0.41±0.11 


1.66±0.18 


1.923 


0.40±0.10 


1.64±0.16 


1.923 


0.35±0.10 


1.56±0.18 


1.923 


0.34±0.10 


1.54±0.18 


1.923 


0.34±0.10 


1.54±0.18 


1.923 


0.37±0.10 


1.59±0.18 


1.82 


0.39±0.10 


1.62±0.18 


2.02 


0.39±0.10 


1.62±0.18 


no tail 


0.38±0.07 


l.eiibO.ll 
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TABLE V. The systematic uncertainties in extracting the ^-u correlation coefficient from 
the Dopplcr Uncshapc of the 35ci(3/2+ , 0) ''^''8(3/2+ , 2939) transition. Note we used a 
Gaussian instrumental lineshape with a centroid E'o, width a, and small low-energy tail. Different 
polynomials of different powers were used to explore the sensitivity to the parameterization of the 
background {i.e. Y^^x^). See text for details. 



parameter amount correlation a 

varied changed 





i = 0.0 


-0.44+0.04 




i = 1.0 


-0.44+0.04 




i = 2.0 


-0.44+0.04 


a 


+20.0% 


-0.31+0.04 


a 


+10.0% 


-0.37+0.04 


a 


-10.0% 


-0.49+0.03 


a 


-20.0% 


-0.55+0.03 


Eo 


+0.4 ch. 


-0.33+0.04 


Eo 


+0.2 ch. 


-0.36+0.04 


Eo 


-0.2 ch. 


-0.49+0.04 


Eo 


-0.4 ch. 


-0.54+0.04 


tail 


+50.0% 


-0.39+0.03 


tail 


-50.0% 


-0.42+0.03 
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TABLE VI. The systematic uncertainties in extracting the ^-u correlation coefficient from 
the Dopplcr lincshapc of the ^501(3/2+ , 0) ''^■''8(5/2+ , 3421) transition. Note wc used a 
Gaussian instrumental lineshape with a centroid Eq-, width cj, and small low-energy tail. Different 
polynomials of different powers were used to explore the sensitivity to the parameterization of the 
background {i.e. Y^^x^). See text for details. 



parameter amount correlation a 

varied changed 





i = 0.0 


-0.39+0.04 




i = 1.0 


-0.38+0.04 




i = 2.0 


-0.39+0.04 


a 


+20.0% 


-0.41+0.04 


a 


+10.0% 


-0.40+0.04 


a 


-10.0% 


-0.39+0.03 


a 


-20.0% 


-0.39+0.03 


Eo 


+0.4 ch. 


-0.39+0.04 


Eo 


+0.2 ch. 


-0.39+0.04 


Eo 


-0.2 ch. 


-0.38+0.04 


Eo 


-0.4 ch. 


-0.39+0.04 



TABLE VIL Toy model results for the ^-v correlation coefficient a in the ^^Cl(3/2^,0) 



^^S(5/2+,3421) 35s(3/2+,0) sequence. The results labeled 'full' corresponds to the complete 
model calculation and the results labeled 'approx' corresponds to a model calculation omitting the 
second forbidden matrix elements. The quoted values are for = —1.26. 



9p 


a 


a 




'fuh' 


'approx' 





-0.240 


-0.216 


-4 


-0.236 


-0.210 


-8 


-0.232 


-0.203 


-12 


-0.228 


-0.196 
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FIG. 1. ^^S energy level diagram showing cleanly identified 7-rays. The energies, spin-parities, 
lifetimes and branching ratios were taken from Ref. [20]. 
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FIG. 2. Time spectra for the 1991 7-ray (bottom), 2939 7-ray (center), and 3421 7-ray (top). 
The points are the experimental data and the sohd hnes are the least squares fits. The muon 
disappearance rate was "divided out" to better demonstrate the hyperfine effect. Note the very 
large hyperfine effect for the 1991 keV 7-ray was used to fix the hyperfine transition rate A/j (see 



Sec. fll F| for details). 
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FIG. 3. Sample fit to the Doppler broadened energy spectrum of tlie 2939 keV 7-ray (tlie 
narrower Gaussian-sliaped peak is a '^^S background line at 2935 keV). 
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FIG. 4. Sample fit to the Doppler broadened energy spectrum of tlie 3421 keV 7-ray (tlie 
weaker Doppler broadened peak is a ^^Ne background line at 3432 keV). 
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FIG. 5. The 7-1^ angular correlation coefficient a versus tlie ratio rjts for the ^^S 2939 keV 
7-ray. Only the statistical errors in a are plotted. 
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FIG. 6. The 7-1^ angular correlation coefficient a versus tlie ratio rjts for the ^^S 3421 keV 
7-ray. Only the statistical errors in a are plotted. 
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FIG. 7. The calculated distribution of ^^Cl ^ ^^S allowed GT strength distribution (BgT) 
using the Is-Od shell model and the universal SD interaction (see text for details). 
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FIG. 8. The calculated capture rate A = 0.12A+ + 0.88A_ for the ^^Cl(3/2+,0) 
3^8(3/2+, 2939) transition (top) and the 35cl(3/2+,0) ^^8(5/2+, 3421) transition (bottom). The 
solid lines are for Qa = —1.26 and the dashed lines are for ga = —1.00. The shaded bands correspond 
to our experimental results. 
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FIG. 9. The hyperfine dependence A+/A- for the =^^01(3/2+, 0) ^ 253(3/2+^2939) transi- 
tion (top) and the 2^Cl(3/2+,0) ^^8(5/2+, 3421) transition (bottom). The sohd hnes are for 
ga = —1.26 and the dashed hnes are for ga = —1.00. The shaded bands correspond to our experi- 
mental results. 
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